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ABSTRACT 



Context. We observe a solar jet at north polar coronal hole (NPCH) using SDO AIA 304 A image data on 3 August 
2010. The jet rises obliquely above the solar limb and then retraces its propagation path to fall back. 
Aims. We numerically model this observed solar jet by implementing a realistic (VAL-C) model of solar temperature. 
Methods. We solve two-dimensional ideal magnetohydrodynamic equations numerically to simulate the observed solar 
jet. We consider a localized velocity pulse that is essentially parallel to the background magnetic field lines and initially 
launched at the top of the solar photosphere. The pulse steepens into a shock at higher altitudes, which triggers plasma 
perturbations that exhibit the observed features of the jet. The typical direction of the pulse also clearly exhibits the 
leading front of the moving jet. 

Results. Our numerical simulations reveal that a large amplitude initial velocity pulse launched at the top of the solar 
photosphere produces in general the observed properties of the jet, e.g., upward and backward average velocities, height, 
width, life-time, and ballistic nature. 

Conclusions. The close matching between the jet observations and numerical simulations provides first strong evidence 
for the formation of this jet by a single velocity pulse. The strong velocity pulse is most likely generated by the low- 
atmospheric reconnection in the polar region which results in triggering of the jet. The downfiowing material of the 
jet most likely vanishes in the next upcoming velocity pulses from lower solar atmosphere, and therefore distinctly 
launched a single jet upward in the solar atmosphere is observed. 

Key words. Magnetohydrodynamics (MHD) - Sun: atmosphere- Sun : corona 



1. Introduction 

[ Solar jets arc the typical collimated and short lived tran- 
. sient ejecta in the solar atmosphere that are significant in 
' mass and energy transport at short spatio-temporal scales. 
, Based on the spatial sizes, life times, and transient nature, 
• the jets arc classified in the form of various phenomena, e.g., 

■ anemone jets (Shibata ct al. 2007), EUV jets (Inncs et al. 
, 1997), spicules (Bohlin et al. 1975), X-ray jets (Shibata et 

■ al. 1992), and penumbral jets (Katsukawa et al. 2007). The 
[ solar surges are cool jets typically formed with a plasma 
. usually evident in and other chromospheric lines, and 
' they are supposed to be triggered by magnetic reconnec- 
I tion processes (Brooks et al. 2007). However, solar surges 

are mostly associated with the flaring regions and the places 
of transient activities at the Sun. Georgakilas et al. (1999) 
have also observed the polar surges which were the cool 
surge-like jets over polar caps formed without any tran- 
sient phenomenon. However, they were having association 
with macro-spicules, and due to the short-lived nature they 
could not attain the greater heights in off-limb polar corona. 

Other short spatio-temporal scale jets (e.g., dynamic 
fibrils, mottles, and spicules) are also significant low- 
atmospheric phenomena to understand its energy and mass 
transport. Using high-resolution observations from 1 m 
Swedish Solar Telescope and advanced numerical model- 
ing, Hansteen et al. (2006) have shown that in the active 



regions these jets are a natural consequence of upwardly 
propagating slow-mode magnetoacoustic shocks generated 
due to the leakage of p-mode sub-photospheric oscillations. 
The active region fibrils are also most likely formed by chro- 
mospheric shock waves that occur when convective flows 
and global oscillations leak into the chromosphere along 
the field lines of magnetic flux concentrations (De Pontieu 
et al. 20G7a). De Pontieu et al. (2004) have also reported for 
the first time the quasi-periodic rise and fall of the spicule 
material due to the leakage of p-mode oscillations from the 
sub-photospheric regions. De Pontieu et al. (2007b) have 
pointed out that typical type-I spicules are found to be 
driven by shock waves that form when global oscillations 
and convective fiows leak into the upper atmosphere along 
magnetic field lines on 3-7 minute timescales. They have 
also found the signature of typc-II spicules, which arc much 
more dynamically formed within 10 s and arc very thin 
(200 km wide) with the lifetimes of 10-150 s. These type 
of spicules seem to be rapidly heated to transition region 
temperatures and eject material through the chromosphere 
at speeds of 50-150 km s~^ . The properties of Type II 
spicules suggest a formation process that is a consequence of 
magnetic reconnection, typically in the vicinity of magnetic 
fiux concentrations in plage and network regions. Recently, 
Morton et al. (2011) have reported a unique signature of the 
relationship between spicule-like structures and the polar 
jet formation using the SDO/AIA observations. However, 
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Fig. 1. The field-of-view of the north polar coronal hole (NPCH) as observed by SDO/AIA 304 A on 03 August 2010 at 02:53:51 
UT. The (Xcen, Ycen) was (-29.618", 918.631"), while the full FOV was (1094', 317") wide with a resolution of 0.6" per pixel. The 
black box shows the triggering of a solar jet that exhibits an exceptional phenomena of its most probable formation by a velocity 
pulse. 



the relationships between various types of jets are still not 
well explored. 

In various types of these jet phenomena, the coronal 
and polar coronal jets are very classic transient phenom- 
ena that have been observed from Skylab and Yohkoh 
era in soft X-rays (Shimojo & Shibata 2000), followed by 
Solar and Heliosphcric Observatory (SoHO) (Wang et al. 
1998), Transition Region and Coronal Explorer (TRACE) 
(Alexander & Fletcher 1999), and Hinode (Kamio et al. 
2007 ; Filippov et al. 2009). The magnetic reconnection is 
found to be one of the most likely mechanisms for the gener- 
ation of typical coronal jets (Shibata ct al. 2007), while the 
reconnection generated waves may also be present with the 
jets during their propagation (Yokoyama & Shibata 1995, 
1996). Pariat et al. (2009) have recently simulated the prop- 
agation of kink motion along the magnetic field lines around 
a rotating and closed magnetic field region. This dynamics 
is most likely associated with the kink instability in the 
jet formation region, which may lead its upward motion in 
the model solar atmosphere. The transversal waves have 
also been recently detected in association with the solar jet 
events at small spatial scales (e.g.. He et al. 2009; Kamio et 
al., 2010, and references cited there). However, the evidence 
of Alfven waves has also been detected in the large-scale po- 
lar X-ray jets (Cirtain et al., 2007). Therefore, the magnetic 
reconnection and associated waves may be the most likely 
drivers to energize such type of the jets at various spatio- 
temporal scales. However, the exact mechanism is not well 
established in the context of the very complex motion of 
the various types of solar jets. 

Recently, Murawski & Zaqarashvili (2010) have re- 
ported the formation of solar spicules by a pulse launched 
within the chromosphere. Their numerical simulations show 
that the strong initial pulse may lead to the quasi-periodic 
rising of the chromospheric material into the lower corona 
in the form of spicules. They have explained the observed 
speed, width, and heights of type I spicules, as well as 
observed multi-structural and bi-directional flows with a 
quasi-periodic rise and fall at 3-5 min time scales due to 
the consecutive shocks. The reconnection generated veloc- 



ity pulses may also drive the solar jets. However, such mod- 
els along with their observational supports are not yet avail- 
able in the context of solar jets. Therefore, it may be worth 
to explore the launching of solar jet due to the formation 
of a pulse in the lower solar atmosphere. 

In the present paper, we find a unique observational sig- 
nature of a jet in the north polar coronal hole which shows 
its dynamics as driven by a velocity pulse. We numerically 
simulate the similar kind of jet triggered by a strong am- 
plitude velocity pulse in the lower solar atmosphere, which 
in general closely mimics the observed properties of the jet. 
In Sec 2, we describe the observational data and analyses. 
We report the numerical model in Sec 3, and results of nu- 
merical simulation in Sec 4. We present the discussion and 
conclusions in the last section. 



2. Observations and results 

We use a time-series data of a solar jet at north polar 
coronal hole as observed in 304 A filter of Atmospheric 
Assembly Imager (AIA) onboard the Solar Dynamics 
Observatory (SDO) on 3 August 2010 during 02:41:14 UT- 
03:02:50 UT. The SDO/AIA has a typical resolution of 0.6" 
per pixel and the highest cadence of 12 s, and it observes 
the full solar disk in three UV (1600 A, 1700 A, 4500 A) and 
seven EUV (171 A, 193 A, 211 A, 94 A, 304 A, 335 A, 131 
A) wavelengths (Lemen et al. 2011). Therefore, it provides 
the unprecedented observations of multi-temperature, high- 
resolution, and high-temporal plasma dynamics all over the 
Sun. The field-of-view of the north polar coronal hole as ob- 
served by SDO/AIA 304 A on 03 August 2010 at 02:53:51 
UT is shown in Fig. 1. The (Xce„, Yce„) was (-29.618", 
918.631"), while the full FOV was (1094", 317") wide with 
the resolution of 0.6" per pixel. The black box shows the 
triggering of a solar jet off the polar limb. The time series 
has been obtained by the SSW cutout service at LMSAL, 
USA, which is corrected for the flat-field and spikes. We run 
aia_prep subroutine of SSW IDL also for further calibration 
and cleaning of the time series data. 
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Fig. 2. The SDO/AIA 171 A (left) and 211 A (right) coronal images overlaid by SDO/AIA 304 A intensity contours. It is clearly 
evident that the coronal filters does not reveal any sensitivity to the jet plasma, which is observed in the He II 304 A filter. 



We observe a very unique jet formation in the east-north 
limb of the polar coronal hole, which seems to be driven by 
a velocity pulse launched in the lower solar atmosphere. 
The jet has been launched slightly obliquely from the back- 
ground open field lines of the polar coronal hole (cf.. Fig. 
2). The life-time of the jet was observed as ^--^21 min, which 
is the typical life time of the coronal jets previously ob- 
served in the Sun and cataloged by Nistico et al. (2009). 
We observe a typical jet that reaches in the corona off the 
polar limb. However, the jet material was only visible in 
He H 304 A filter, which is sensitive to the comparatively 
cool temperature of 1x10^ K. The jet material could not 
be evident in the typical TR/coronal filters (e.g., Fe IX 171 
A , 211 A filters) of SDO/AIA, which are sensitive to the 
plasma temperature >10^ K. Fig. 2 shows the SDO/AIA 
171 A (left) and 211 A (right) coronal images overlaid by 
SDO/AIA 304 A intensity contours at a time of the max- 
imum rise of the jet. It is clearly evident that the coronal 
filters do not exhibit any sensitivity to the jet plasma, which 
is only observed in the He II 304 A filter. This scenario per- 
sist during whole life-time of the observed jet. The unique 
jet is a very faint event even captured in 304 A as most 
of the strong emissions come from spicule borders at this 
wavelength also. However, we may not rule out some very 
weak amount of coronal emissions also that could not be 
distinguished in the diffused background off the polar limb. 
Therefore, we abbreviate this jet as a cool jet as its plasma 
mostly emits the radiations captured by the AIA filters 
that are sensitive to cooler plasma emissions, e.g., 304 A. 
There may also be the possibility of some traces of emis- 
sions m other SDO/AIA filters (e.g., 1600 A ) that are sen- 
sitive to the emissions from comparatively cooler plasma. 
This means that the jet material was typically formed by 
comparatively cooler plasma compared to the background 
coronal plasma, and it moves in form of plasma column 
that does not exchange any heat with the ambient coronal 
plasma. This was the unique property of the observed cool 
jet off-limb in the polar coronal hole. Our observed cool 
ejecta reaches at a higher height (72 Mm) with a life-time 



of 21 min and larger width of 20 Mm. We, therefore, ab- 
breviate it as a cool polar jet, rather as an evidence of any 
polar surge like ejecta. 

The jet starts moving up above the solar limb on ^02:41 
UT and reaches at a maximum height of ~72 Mm at 02:52 
UT. Therefore, its average rising speed is ^110 km s~^ that 
is the typical speed of the jet by which it rose in the corona 
at a maximum height of ~72 Mm in the first ^^648 s of the 
time (Fig. 3). Although, the width slightly varies during the 
course of time at the base of the jet, its shape is approxi- 
mately uniform during its rise-up. Its width at the base is 
measured approximately as ~20 Mm. The unique evidence 
this jet presented is its rise straight in the corona off-limb 
in first 10-11 min of its life time, and then its fall back from 
the same path in rest of the duration of its life time. Some 
unique observations of chromospheric jet by Hinode/SOT 
Ca II H filter revealed fine structuring as well as transver- 
sal motion that was associated with the helicity injection 
in the upper solar atmosphere (Liu et al. 2009). However, 
the spatial resolution of SDO/AIA (0.6") is indeed lower 
compared to SOT (0.2"), and secondly the observed cool 
jet is very fainter to clearly identify such fine structuring. 
Moreover, the observed intriguing chromospheric jet (Liu 
et al. 2009) was associated with a small flare at the western 
limb, while in the present case we observe the cool jet in 
polar coronal hole. Therefore, these two cases are entirely 
different, and we do not observe such type of oscillating 
threads as well as helicity injection in our observations. 

The similar type of the dynamics in form of the 
parabolic trajectory has been reported in case of various 
types of spicules (type I and II), dynamic fibrils, mottles 
(Hansteen et al. 2006; De Pontieu et al. 2007a,b, 2011). 
These findings also support the acceleration of such small- 
scale jets by the shock triggered in the chromosphere due to 
velocity pulses generated by the leakage of sub-photospheric 
p- modes. However, such drivers need few km s^^ velocity 
amplitude, 3-5 minutes temporal scales, as well as some con- 
centrated field regions of few hundred Gauss to push such 
jets 2-10 Mm above the solar limb. Recently, Murawski & 
Zaqarashvili (2010) have also simulated the solar spicules 
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Fig. 3. The SDO/AIA 304 A time sequence of a cool jet propagation above the solar limb in the north-east polar coronal hole 
that clearly shows the upward motion of the plasma and its downward motion along the same path. 
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by launching a velocity pulse above the solar photosphere. 
However, the strong velocity pulse will be needed in the 
present case to trigger a jet material higher in the corona. 

Our observed dynamics supports that this jet may also 
be formed by a velocity pulse that triggered in the chro- 
mosphere and steepens in form of a shock in the transition 
region/corona. However, it may not be associated with typ- 
ical photospheric velocity pulses of the order of few km s^^ 
associated with the leakage of p-modes that launches the 
small-scale jets as mentioned in many significant previous 
studies (Hansteen et al. 2006; De Pontieu et al. 2007a, b). 
In the present case, the steepened shock can be generated 
by a reconnection generated strong velocity pulse. Its front 
may rarefy the plasma density behind it, and the jet plasma 
could move up. After reaching in the corona along a curved 
magnetic field lines, this jet falls back due to gravity and 
cool plasma traced back its backward path as observed here. 
In next section, wc discuss in details the numerical simula- 
tion of the observed jet in the solar atmosphere triggerred 
by a strong velocity pulse. 



3. A numerical model 

Our model system is taken to be composed of a 
gravitationally-stratified solar atmosphere which is de- 
scribed by the ideal two-dimensional (2D) MHD equations: 





9q „ 
— + V 




= 0, 


(1) 


V)V = 


-Vp+ -(V X B) 
M 


X B + 


eg, 


(2) 


dp 

at 
In ' 


+ V ■ (pV) = (1 - 


- 7)pV 


■ V, 


(3) 


= V X (V X B) , 


V B 


0. 


(4) 



Here g is mass density, V is flow velocity, B is the magnetic 
field, p = ^qT is gas pressure, T is temperature, 7 = 5/3 
is the adiabatic index, g = (0, —g) is gravitational acceler- 
ation of its value g = 274 m s"'^, m is mean particle mass 
and fee is the Boltzmann's constant. 

3.1. Equilibrium configuration 

We assume that at its equilibrium the solar atmosphere is 
still (Vc = 0) with a force-free magnetic field. 



(V X Be) X Be = 



(5) 



such that it satisfies a current-free condition, V x ~ 0, 
and it is specified by the magnetic flux function, as 

Be = V X (Az) . 

Here the subscript e corresponds to equilibrium quantities. 
We set an arcade magnetic field by choosing 



A{x, y) = BqAb cos (x/AB)exp[-(y - j/iO/Ae 



(6) 



Here. Bq is the magnetic field sX y ~ yi-, and the magnetic 
scale-height is 



Ab = 2L/tt . 

We set and hold fixed L = 100 Mm. 



(7) 



As a result of Eq. (5) the pressure gradient is balanced 
by the gravity force, 



Vpe + Peg = . 



(8) 



With a use of the ideal gas law and the i/-component of Eq. 
(8), we arrive at 



Pc[y) = Po exp 
where 



dy 

'v. My') 



Pcjy) 
gHy) ' 



A{y) = fcBTe(y)/(m.g) 



(9) 



(10) 



is the pressure scale-height, and po denotes the gas pressure 
at the reference level that we choose in the solar corona at 
y, = 10 Mm. 

We adopt an equilibrium temperature profile To(z) for 
the solar atmosphere that is close to the VAL-C atmo- 
spheric model of Vernazza et al. (1981). Then with the 
use of Eq. (9)we obtain the corresponding gas pressure and 
mass density profiles. 

3.2. Perturbations 

We initially perturb the above equilibrium impulsively by 
a Gaussian pulse in the vertical component of velocity, Vy, 
viz., 



Vy{x, y,t = 0) = Av exp 



(x - xo)"^ (y - 2/o)^ 



(11) 



Here A^ is the amplitude of the pulse, (xo,yo) is its initial 
position and lUx, Wy denote its widths along the x- and y- 
directions, respectively. We set and hold fixed A^ = 160 km 
s-\ xo = -17.5 Mm, yo = 1.75 Mm, = 10 Mm, and 



Wy = 0.1 Mm. 



4. Results of numerical simulations 

Equations (l)-(4)are solved numerically using the code 
FLASH (Lee & Deane 2009). This code implements a 
second-order unsplit Godunov solver with various slope lim- 
iters and Riemann solvers, as well as adaptive mesh refine- 
ment (AMR). We set the simulation box of (—125, 15) Mmx 
(0, 110) Mm along the x- and y-directions and impose fixed 
in time all plasma quantities at all four boundaries of the 
simulation region. In all our studies we use AMR grid with 
a minimum (maximum) level of refinement set to 4 (8). 
The refinement strategy is based on controlling numerical 
errors in mass density, which results in an excellent resolu- 
tion of steep spatial profiles and greatly reduces numerical 
diffusion at these locations. 

Figure 4 displays the spatial profiles of plasma tempera- 
ture (colour maps) and velocity (arrows) resulting from the 
initial velocity pulse of Eq. (11), which splits into counter- 
propagating parts. The part which propagates downwards 
becomes reflected from the dense plasma layers at the pho- 
tospheric region. This reflected part lags behind the origi- 
nally upward propagating signal which becomes a shock. As 
the plasma is initially pushed upwards the under-pressure 
results in the region below the initial pulse. This under- 
pressure sucks up comparatively cold chromospheric plasma 
which lags behind the shock front at higher temperature. 
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Fig. 4. Temperature (colour maps) and velocity (arrows) profiles at t = 25 s, t = 200 s, t = 500 s, t = 600 s, t = 700 s, t — 800 
s, t = 900 s, t — 1000 s, and t = 1200 s (from top-left to bottom-right). Temperature is drawn in units of 1 MK. The arrow below 
each panel represents the length of the velocity vector, expressed in units of 150 km s^^. 
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As a result, the pressure gradient force works against grav- 
ity and forces the chromospheric material to penetrate into 
the solar corona in form of rarefaction wave. At i = 25 s 
this shock reaches the altitude of ?/ ~ 10 Mm and the rar- 
efaction wave is located 1 Mm below the slow shock. The 
next snapshot (top middle panel) is drawn for t = 200 s. 
At this time the slow shock reached the altitude of y ~ 63 
Mm while the chromospheric plasma is located below at 
y ~ 40 Mm. At the next moment of time, t — 500 s, (top 
right panel) the shock moved up and the chromospheric 
plasma blob exhibits its developed phase, reaching the level 
of 2/ = 64 Mm, which matches well with the observational 
data of Fig. 3(top right panel). The cool jet slow down while 
propagating upwards. At t ~ 800 s it arrived at the level 
of y ~ 70 Mm (Fig. 3, middle-right panel)and subsequently 
subsided as the plasma began to flow downward, being at- 
tracted by gravity (bottom panels). At the subsiding stage 
the fine structures developed as a result of interaction be- 
tween falling down plasma and rising up secondary pulse. 
Similar features were evident in the case of the numerical 
simulation of spicules (Murawski & Zaqarashvili 2010). 

5. Discussions and Conclusion 

Our numerical simulations show a general scenario of ris- 
ing and subsequently falling down cool plasma jet with a 
number of observational features like jet height and width, 
ballistic motion, approximate rising time-scale and aver- 
age rising velocity that are being close to the observational 
data. However, some match between the numerical and ob- 
servational data is approximate and rather qualitative, e.g., 
the rising time-scale in the simulation is ^800 s while in the 
real observations it is ~720 s. This mismatch may obviously 
result from some of the simplified profiles of the model pa- 
rameters in our developed numerical model, e.g., the mag- 
netic field configuration at a place at which the initial pulse 
was launched, and so on. Moreover, the real jet was excited 
in more complex plasma and magnetic field conditions at 
the polar coronal hole of the Sun. We might not model ex- 
actly realistically the way in which the jet was excited in 
the real Sun. In fact, the exciter could work for some time, 
it could be located at a different place and it could have a 
different size, speed and distribution. As such data is not 
provided by the observations we decided to trigger the jet 
by a localized pulse in a plasma velocity, that is launched 
below the transition region. By this way we managed to ex- 
cite the jet which mimics on average the properties of the 
observed cool plasma jet. The small mismatches like the 
exact height and length of the jet as well as its descend- 
ing time can be reduced by tuning of free parameters in 
the numerical model. In particular, the jet length would be 
smaller for a smaller amplitude of the initial pulse. We did 
not expect to fit perfectly our numerical data to the ob- 
servational findings but instead our intention was to show 
a qualitative and to some extent quantitative agreement 
between the numerical and observational approaches. 

It is noteworthy that recently Culhane ct al. (2007) have 
found the first observational evidence of unique EUV polar 
jets which could not escape from the Sun as like the jets 
previously observed with SXT in Yohkoh era. They con- 
ject that such accelerated polar jet plasma is being heated 
upto coronal temperature with its rise, and then fall back to 
lower atmosphere after its cooling. They observed jet veloc- 
ities in the range of 150 km s~^- 360 km s~^ and reported 



their formation most probably due to the magnetic recon- 
nection between emerging bipoles with the large-scale open 
field lines of the polar coronal hole. Ko et al. (2005) have 
also found cool and hot plasma components in a limb jet, 
as well as its falling back scenario in the lower atmosphere. 
However, all these previously observed jets were found to 
be consistent with the impulsive transient heating (e.g., 
as driven by direct reconncction processes), and thereafter 
cooling and draining along the same path. We also observe 
the jet motion that moves up and thereafter falls along the 
same path in the lower atmosphere. However, the jet is 
made by the cool plasma maintained at a temperature of 
--10^ K, which is sensitive to the SDO/AIA 304 A fiher rep- 
resenting upper chromospheric/TR plasma. However, the 
jet plasma column is not evident at typical coronal temper- 
ature of 1.0 MK in the SDO/AIA 171 A . Therefore, the jet 
consists of a cool plasma even after reaching at the coronal 
heights, and does not subject to any heating and cooling 
events as previously observed by Culhane et al. (2007); Ko 
et al. (2005). The primary mechanism for the origin of such 
jet may be the magnetic reconncction between the emerging 
bipolar magnetic fields and pre-existing open field lines in 
the polar coronal hole (Culhane et al. 2007). However, com- 
paratively slow average rising speed (110 km s~^), steady 
ballistic motion, the absence of impulsive heating at the 
base and no gradual enhancement in plasma temperature 
etc., indicate the role of magnetic reconncction as an in- 
direct mechanism. The jet is seem to be driven by a large 
amplitude velocity pulse that may be excited non-linearly 
by recurrent magnetic reconncction at the base of the po- 
lar coronal hole. This pulse steepens into a slow shock front 
which moves at a higher temperature wave front. However, 
the density and typical emissions are low in polar coronal 
holes. Therefore, we could not observe high temperature 
bright front of the jet. Later, the motion of the slow shock, 
under stable conditions, causes the propagative rarefaction 
wave which direct the cool and denser chromospheric/TR 
plasma upward. This steadily moving plasma is clearly vis- 
ible in the observations as well as in the simulation. After 
reaching at the height of ^75 Mm, the plasma begins to 
subside. The gravitational free fall plays the role for the 
falling of jet material. During the subsiding phase of the 
jet, in the numerical simulation, the another pulse is evi- 
dent at the base that causes the low rate of subsided plasma 
compared to the observations. However, the realistic solar 
atmosphere is rather complex, and there is a possibility 
that the huge downfalling material causes the vanishing of 
the another pulses coming from the lower atmosphere. This 
is the reason why a pulse driven single jet became evident 
above the limb rather than multiple jets at the same place. 

In conclusion we suggest that the initial velocity pulse 
launched below the transition region is able to trigger a 
hot plasma shock which is followed by a cool plasma jet. 
This cold plasma jet approximately resembles many of the 
features which are found in the presented SDO/AIA ob- 
servational data. We report first time on the observations 
of a pulse driven plasma jet in the polar region and pro- 
vide theoretical explanation of this phenomenon on the ba- 
sis of numerical simulations we performed. However, fur- 
ther multiwavelength observations should be performed by 
high-resolution space borne (e.g., SDO, Hinode, STEREO) 
and complementary ground based observations to shed new 
light on this kind of unique jets that rise and fall back in 
the solar atmosphere. This will also impose a rigid con- 
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straint on the stringent simulations of such kind of jets in 
the model solar atmosphere. 
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